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Abstract 

In this work, masses and radii of neutron stars are considered to investigate the effect of nuclear 
symmetry energy to the astrophysical observables. A relativistic mean field model with density- 
dependent meson-baryon coupling constants is employed in describing the equation of state of dense 
nuclear matter, and the density dependencies of the symmetry energies are quoted from the recent 
phenomenological formulae obtained from the heavy ion data at subnuclear saturation densities. 
Since hyperons can take part in the /3-equilibrium of the dense matter inside neutron stars, we 
include hyperons in our estimation and their roles are discussed in combination with that of the 
nuclear symmetry energy. 

PACS numbers: 13.75.Ev, 21.65.Ef, 26.60.Kp 
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I. INTRODUCTION 



Recent works about the nuclear symmetry energies at densities around the nuclear sat- 
uration reveal an interesting diversity 7t~/tt + production ratio in FOPI data of the heavy 
ion collisions at SIS/GSI favors super-soft behavior of the symmetry energy, which even 
becomes negative at densities above a few times the saturation density [l). On the other 
hand, analysis of the giant monopole resonance for various nuclei, and the determination 
of the neutron skin thickness show a monotonic increase of the nuclear symmetry energy 
as density increases, even though the stiffness among them is substantially different (2), yj. 
Current status seems to demand more theoretical analysis as well as ample experimental 
data in order to sharpen our understanding on the nature of nuclear symmetry energy. 

Neutron stars can provide a benchmark to test the density dependence of the symmetry 
energy. Analysis for the neutron star with various forms of symmetry energy was performed 
with a simple free gas model Q , and the result for the matter composition in the interior of 
the neutron star shows high sensitivity to the density dependence of the symmetry energy. 
A softer symmetry energy allows more neutrons in the neutron star matter. On the other 
hand a stiffer one requires more energy for the matter with higher asymmetry, and thus it 
favors larger proton fraction. With a very soft equation of state (EOS), it is difficult for the 
neutron star matter to sustain the strong gravitational contraction, and the star will easily 
collapse to a black hole. Conversely, a stiffer EOS allows a larger value of the maximum 
mass of a neutron star. Well-measured neutron star masses lie in the range of (1.2 ~ 2.0)M Q 



6] and both super-soft and super-hard EOS which give neutron star masses far outside 
of this range can be excluded . However, the maximum mass of neutron star is still an open 
question. 

In this work, our main interest is the role of symmetry energies which are fixed from the 
analysis of modern experiments to the bulk properties of neutron stars. It has been shown 
in numerous works that hyperons play a crucial role in making the neutron star matter 
significantly soft and consequently reduce the maximum mass of the neutron star [7 13]. 



When only nucleons are considered, soft symmetry energy will produce a matter exclusively 
dominated by the neutron, but if the flavor changing /3-equilibrium is allowed, sufficiently 
high Fermi moment of the neutron makes the weak decay such as n — > A much more feasible. 
These decays will have a critical impact on the composition, cooling rate, and EOS of the 
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neutron star matter. Therefore it is essential to include hyperons in the neutron star matter. 

Production of hyperons and/or transition to exotic states such as quark deconfmement or 
Bose-Einstein condensation are known to be very sensitive to the input parameters such as 
coupling constants and their density dependencies. In Ref. [ijj], we considered the density- 
dependent masses and coupling constants for u and p mesons, and showed that the density 
dependence affects the equation of state and bulk properties of neutron stars critically. In 
this work, we employ four empirical formulae for the symmetry energy and plug them into a 
mean field model in which the coupling constants are density-dependent. By adopting such 
a hybrid method, we can single out the effect of the symmetry energy on top of well-defined 
nuclear saturation properties, and investigate its role to observables. 

The empirical formulae for the symmetry energy are determined from the data at sub- 
nuclear saturation densities in the relativistic heavy ion experiments. Direct information 
about the behavior of the symmetry energy at densities above the saturation is still lacking. 
However, by extrapolating the existing formulae to supranuclear densities, and calculating 
the neutron star properties with them, we can indirectly constrain the behavior of the sym- 
metry energy at high-density region. The results indeed show that the different behavior 
gives measurable differences in the properties of neutron stars. 

The paper is organized as follows. In Sec. II, we introduce the models and formulae used 
in the calculation. In Sec. Ill, numerical results are displayed and we discuss them. We 
conclude the paper in Sec. IV. 



II. MODEL AND SYMMETRY ENERGY 



A. Density-dependent relativistic mean-field model 



We briefly introduce the density-dependent relativistic mean-field (DDRMF) model with 



octet baryons, and a, u>, and p mesons 
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-|l9j. The model Lagrangian is given as 



C = J> 



T 



(m B - T aB a) 



if, 



B 



+ l - {d^a - mla 2 ) - + ^uy^ - ■ IT + ±m% ■ fr, (1) 

where I\b are density- dependent meson-baryon coupling constants, and u' IU (= d^u u — d u u fJ- ) 
and b tJlU (= d^b u — d v U L ) are the field tensors for uj and p meson fields, respectively. 



For the infinite nuclear matter, the equations of meson fields in the mean-field approxi- 
mation are obtained as 



B 



bPb 



<bPb 



B 



m 



^03 — ^ F pB I m p B , 



(2) 



where p B (= ip B ip B ) and ps(= i^ b i^b) ar e the scalar and vector densities of a baryon £>, 
respectively, and I B3 is its isospin ^-component. The Dirac equation of a baryon can be 
written as 



(3) 



If we assume T iB = T iB (p) where p = ^2 b Pb, the scalar and time component of vector 
self-energies are obtained as 



e* = r ff Bff 

s° = r„ B u + r pB b 03 i m + z v 



with 



E 0(r) = 



dT^B . ^r^B r , dT aB _, 

— — WoPb + -5 — Ib3Pbo 3 n — o-ps 
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(4) 



(5) 



Note that we only need time component of in the mean field approximation. 

With the above equations, one can obtain the energy density (e) and the pressure (P) in 
the form 

£ = YsJ^y // d 3 ke B (k) + -[mla 2 + mlu 2 +m 2 p b 2 03 \ 



b 



(6) 



where, eg(A;) = ^Jk 2 + (m B ) 2 and is the effective mass of a baryon given by m B = 
m B — T aB a. Note that the rearrangement self-energy term £? contributes to the pressure 
explicitly, but doesn't to the energy density 19j. Chemical potential of a baryon at its Fermi 
surface reads 

de 



Pb 



dpB 



e B (k$) + T^coo + I m T pB b 03 + £°< r >. 



(7) 
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Meson(i) mj(MeV) a, 6j 



a 550 10.87854 1.365469 0.226061 0.409704 0.901995 

uj 783 13.29015 1.402488 0.172577 0.344293 0.983955 

TABLE I: Parameters of the density-dependent coupling constants in Typel & Wolter |.1.7] lil- 
ted to the saturation density po = 0.153 fm -3 , binding energy per nucleon 16.247 MeV, and the 
compression modulus Kq = 240 MeV. 

B. Coupling constants for symmetric matter 

In this work, we employ the density-dependent coupling constants for the a and u mesons 



proposed by Typel & Wolter (TW99) [17|. Coupling constants for the p meson will be 



discussed in the next subsection. For i = a, u>, we can write 

r«£j(p) = gwfi(n), (8) 

where n = p/p with p the saturation density. It is assumed that fi(l) = 1, so gm denotes 
the coupling constant at the saturation density. Density-dependent part fi(n) is given by 

1 + b i (n + d i ) 2 

m = a * i+cAn +dl r (9) 



Details for the procedure to fix the parameters in fi(n) can be found in [13], and we simply 
display the values of the parameters in Table [J and the resulting saturation properties in 
the caption. 

For the couplings between hyperons and cu-meson, we assume the quark counting rule, 
he. g u A = 3wS = §<7u>iV f° r A and S, and g^ = \g^N for S. For the hyperon-cr meson 
coupling constants, we fix them to reproduce the optical potential values for A, £ and H at 
the saturation density, —30 MeV, 30 MeV and —15 MeV, respectively, which are determined 
empirically. Numerical results thus determined are g^y/gaN = 0.627, 0.480, 0.313 for Y = 
A, S and 5, respectively. As for the density dependence of the T^y, we assume the same 
dependence as given by T^, Eq. dHJ) for simplicity. The density dependence and coupling 
constants of p meson will be discussed in the next subsection. 
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C. Symmetry energy beyond nuclear matter density 



Energy per baryon in infinite nuclear matter is conventionally defined as 

E(p, 5) = E(p, 0) + E sym (p) 5 2 + 0(5 4 ), (10) 

where p = p n + p p is the nucleon number density and 5 = (p n — p p ) / p. From the above 
equation, we can define the symmetry energy as 

E sym (p) = E(p, l)-E(p, 0). (11) 

With the DDRMF model described above, the energy density in the nucleonic phase is 
obtained as 

£ = E T^p / ' d 3 ke N (k) + l - [mla 2 + m> 2 + m 2 p b 2 03 \ . (12) 

i=n,p v / 

Since the a and the u fields depend only on the baryon number density p, we obtain 

£ S ym = e(p, 5 = 1)- e(p, 5 = 0) 

If we divide the symmetry energy into kinetic and potential terms as 

-Esym T S y m -\- V^y m , (14) 

the kinetic term reads 

T ^ = {kNpf (15) 

P 6v/(^)2 + (m^' 

where kp is the Fermi momentum of the nucleon in symmetric nuclear matter. The potential 

part is then written as 

v;,, = 3|p. (is) 

In the conventional considerations, T pN is assumed to be density-independent, and its value 
is determined from the empirical symmetry energy at the nuclear saturation density, E sym = 
30 ~ 35 MeV . However, recent works about the symmetry energy show that its density- 
dependence can be much more significant than what has been understood so far. Especially, 
the uncertainties in its high-density behavior range from super-soft to very hard ones. 

Remaining part of this subsection is devoted to the description of the density-dependence 
of the symmetry energies employed in this work. 
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MDI 

Symmetry energy has been considered in various theoretical frameworks such as two- 
and three-body nuclear forces, Brueckner-Hartree-Fock formalism, relativistic mean 



field theories, effective field theories and etc. In Ref. [20[, the authors included 
momentum-dependent interaction (denoted as MDI) in the isovector channel of the 
interactions and obtained the so called MDI version of the symmetry energy. The 
result is parameterized in the form 

E sym (n) = T sym + A{x)n + [18.6 - A{x)]n B ^ x \ (17) 

where x is a parameter fixed from experimental data. Behavior of the symmetry 
energy is sensitive to the parameter x, and three representative values are 1, and 
— 1. Among these values, production ratio tt~/tt + in FOPI data is reproduced well 
with x — 1, with which A(x) ~ 107 MeV and B(x) ~ 1.25 The potential terms 
derived from the above form are related to the one in the mean field model as 

■^?n = A.(x)n + [18.6 - A(x)]n B{x) . (18) 

Since B{x) is different from 1, if we divide both sides with n, n B ^~ 1 remains density- 
dependent, and consequently we have a density dependence for the coupling constant 

FSU & NL3 

Symmetry energy can be expanded in powers of (n — 1). At subnuclear saturation 
densities, symmetry energy can be approximated to quadratic order, and it can be 
written in the form 

E sym (n) = J + -L{n - 1) + — K sym (n - l) 2 , (19) 
o lo 



where J, L and K sjm are the parameters either calculatec 



the measured properties of nuclei. FSU [2l| and NL3 



from theories or fixed from 
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indicate the parameter 



sets of non-linear QHD model. In the NL3 model, cubic and quartic terms of the 
self interactions of a mesons are added to the quadratic QHD model. In the FSU 
model, non-zero contribution of a term proportional to cj^c^ ■ is added to the 
NL3 one. Coupling constants and parameters of the model are fitted to the static 



Model 


Po 


E/A 


K 


J 


L 


K 


FSU 


0.148 


-16.30 


230.0 


32.59 


60.5 


-51.3 


NL3 


0.148 


-16.24 


271.5 


37.29 


118.2 


100.9 



TABLE II: Nuclear satuation density pq is in fm 3 , and the binding energy per nucleon E/A, the 
compression modulus Kq, and the parameters in the symmetry energy J, L and K sjrn are in MeV. 

properties of several nuclei. Once the parameters of the model are determined, one 
can calculate the parameters that enter the approximate form of the symmetry energy 
in Eq. ( Fl9|) . Results for the nuclear saturation properties and the parameters for the 
symmetry energy in Eq. ( TPO!) are shown in Table [Til Data from isospin diffusion and 
giant monopole resonance in heavy ion collisions give the values L = 88 ± 25 MeV 
and K sym — QL = -500 ± 50 MeV. Values of L and K sym from both FSU and NL3 
are within or close to the empirical range. One has to note that the kinetic term is 
not included in Eq. ( fl9l . As a result, the potential term in the FSU and NL3 models 
reads 

T2pNP °- ~ J + l ~L{n - 1) + ±-K sym (n - l) 2 - T sym . (20) 



-n 



8m 2 p " ' 3 v " _/ ' 18 
• TW99 

In Ref . [l7| , the density dependence of T pN is given as 

r p jv = #pjvexp[-a p (n - 1)]. (21) 
A calculation based on Dirac-Brueckner approach showed that the pN coupling con- 



231 ] . The model gives g p ^ 



stant becomes small in magnitude at high densities 
7.32196, a p = 0.515 and the symmetry energy E sym = 33.39 MeV at the saturation 
density. 

In Fig. [TJ we show V sym (p) for the four models. In the MDI model, since A(x) > 18.6 
and B(x) > 1, V^ ym is an increasing function at small n and reaches maximum at n = 0.88. 
After the maximum it decreases monotonically, and becomes negative for n > 2.1. A similar 
trend is observed from the curve for the FSU model. In this model, the symmetry energy 
becomes negative for n > 7.6. 
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FIG. 1: V sym (p) = T 2 pN p/(8m 2 p ) for four models. 

Since we assume V sym (n) = T 2 N p n/8m 2 in our calculation, negative values of V sym imply 
negative values of T 2 pN . In the equations that determine the EOS, T p n always enters in the 
quadratic form, so negative values of T 2 pN cause no mathematical difficulty in solving the 
equations. However, it can be problematic physically because a negative T 2 pN value gives 
a pure imaginary number for the p- meson field, &03- We want to note that in most cases, 
V sym (p) of the models is determined from the data at densities around or less than the 
saturation density. It is not clear yet to what density the extrapolation is valid. At the 
same time, it is not clear either to what density the term r 2 N pon/8m 2 will be dominant in 
the symmetry energy. In other words, not only the quadratic term but also higher order 
terms of the p-meson contribution in various forms can be important at high densities. 
Negative Y 2 pN value in the present work may be regarded as either an unphysical result due 
to naive extrapolation of the empirical formula to high densities, or a result of higher order 
contributions that can give rise to negative values of the symmetry energy. 

In the NL3 model, V^ ym is a quadratic function with a positive coefficient for the quadratic 
term, and thus it is a monotonically increasing function. In the TW99 model, T p n decreases 
exponentially, so V sym converges to zero at high densities. 

At n > 1, the four models show very distinctive dependence on the density, and we may 
classify NL3 to a hard model, FSU and TW99 to a soft one, and the MDI to a super-soft 
one. These controversial behaviors put distinguishing imprints on the properties of neutron 
stars. In next section, we discuss the implications of symmetry energies to the properties of 
nuclear matter inside neutron stars. 
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III. NEUTRON STAR EQUATION OF STATE 
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FIG. 2: Populations of particles for np. 

In this section, we discuss the nuclear matter properties with various forms of symmetry 
energy. For the comparison we considered the nuclear matter with nucleons only (np) and 
with full baryon octet including hyperons (npH). 

A. Chemical composition of nuclear matter 

Figure [2] shows the fractional ratio of particle (nucleon and lepton) densities to the baryon 
density for np. For the MDI model (Fig. [2(a)), the proton fraction is less than 10 % at 
densities where V sym > 0. As V syra decreases quickly at high densities, only the neutrons 
remain in the star. Such an unusual particle composition can be understood from the 
behavior of V^ ym : Negative values of V sym favor more asymmetric matter than the symmetric 
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one. With the proton fraction less than 10 %, the direct URCA process, n — > pe~i? e cannot 
happen in the star. If there are only neutrons, then the weak decay will no longer happen, 
and neutrino emission will cease from the interior of the star. Emission of the neutrino has a 
significant effect to the thermal evolution of the neutron star, and thus, though it is beyond 
the scope of this paper, it is important to investigate whether the symmetry energies are 
compatible with the cooling curve of the neutron star. The particle fraction for FSU model 
(Fig. [2(b)) shows similar qualitative behavior as that of the MDI model in the region where 
the symmetry energies become negative. 

In the TW99 model, V sym converges to zero at high densities, which means that the 
isospin dependent interaction doesn't affect the /3-equilibrium condition any more. Since a- 
and w-mediated interactions are isospin-independent, only the kinetic term determines the 
particle fraction as V^ ym approaches zero in the TW99 model. This is the reason why the 
particle fractions saturate at high densities in this model as one can see in Fig. [2(d). On 
the other hand, V sym increases monotonically with density in the NL3 model (Fig. E(c)). As 
a result the proton fraction increases consistently as density increases. 

In Fig. [31 we show the fractional ratio of particles for npH. One can see that the critical 
densities for the creation of hyperons, and their fractional ratios above the critical densities 
are very sensitive to the density dependence of the symmetry energy. With hyperons, an 
unusual result appears from the MDI model (Fig. 13(a)), which is very contrasting to those 
of the remaining models. A is the lightest hyperon, and at the same time it feels the 
strongest attraction among the hyperons in the baryon octet. For these reasons, A hyperon 
is the first hyperon that are created via /3-equilibrium condition in the FSU, NL3 and TW3 
models (Fig. [3(b), (c), (d), respectively). On the other hand, since the E hyperons feel 
the strongest repulsion among the hyperons in the nuclear matter, their chemical potential 
is always too large to satisfy the /3-equilibrium conditions, so it never appears in the NL3 
model. Similar result was obtained from the quark-meson coupling model 24] • I n the result 
with MDI model, however, even though £ hyperon is repulsive, it is the only hyperon that 
satisfies the /9-equilibrium condition. 

We can understand the unusual behavior in the MDI model from the chemical potentials. 
The chemical potential of n, A and E~ read respectively as, 



fin = e n + T^nujq + - ( ) (pn - p P ), (22) 
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FIG. 3: Populations of particles for npH. 
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In the MDI model, p n — Pp increases in the region where V^m decreases. Even when V sym < 0, 
p n — p p keeps increasing, but at densities where V^ ym is negative, the p-meson contribution 
to p n becomes negative. In other words, in the NL3 model, the p-meson contribution to 
p n is always positive but it becomes negative above a certain density in the MDI model. 
On the other hand, since there is no p-meson contribution to pa, the chemical potential of 
A hyperon increases with density. If the negative contribution of p-meson to p n becomes 
significant before the /3-equilibrium condition p n = pa is satisfied, then there is no chance to 
have A hyperons in the interior of the neutron star. As for the E hyperon, since T p s oc T p n 
in our calculation, when the p-meson contribution to p n is negative, it is also negative to 
Ps-- Since the /3-equilibrium for E~ reads p n + p e - = ps-, if the Fermi momentum of 
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FIG. 4: Equation of state for np (left pannel) and npH (right pannel). 

electrons increases sufficiently to satisfy the /3-equilibrium condition for £~, we can have S~ 
prior to A, as summarized in Fig. |3j 

An early creation of A hyperon in the NL3 model than in the FSU and TW99 models 
is also interesting even though the differences in the critical densities are not so big. The 
/3-equilibrium condition hk = n P + /i e - , and the behavior of symmetry energies provide us 
a simple and useful insight into the origin of the difference. Because V sym in the NL3 model 
is stiffer than those in the FSU and TW99 models, the proton fraction in the NL3 model is 
roughly twice of those in the FSU and TW99 models at densities around n = 2. A larger 
proton fraction gives a larger Fermi momentum of the proton and the electron, hence the 
early creation of A is favored, reducing the energy of the system. 

B. Neutron star equation of state 

In Fig. HI we summarize the pressure as a function of density for both np (Fig. 11(a)) and 
npH (Fig. H(b)). Since we use a single model for the isospin independent interactions, the 
difference in the pressure directly reflects the different contribution from isospin dependent 
interactions, V^ ym as in Eq. <Q. Similar to the curves for V sym in Fig. [TJ the NL3 model 
exhibits the stiffest increase among the models, the FSU and TW99 models show milder 
behavior than the NL3, and the MDI model shows dramatic softening of the EOS due to 
the symmetry energy. In the MDI model, the pressure decreases above p = G.5po, which 
is a signal for the instability. In the results with hyperons, however, the curves are more 
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FIG. 5: The relation between mass and radius for np (left pannel) and npH (right pannel). 

complicated. With a stiffer symmetry energy, we discussed that hyperons are created at 
lower densities. Because the isospin-independent repulsive interaction of hyperons is weaker 
than that of the nucleon, the appearance of hyperons makes the pressure softer in the 
isospin-independent part. Complicatedness of the curves for the pressure with hyperons 
than that with only nucleons may be due to the competition between the enhancement 
from the isospin-dependent interaction and the reduction from the isospin-independent one. 
Comparing the pressure from NL3 to the ones from FSU or TW99, one can see that the 
pressure from NL3 is stiffer at low densities, but it becomes softer at high densities. This 



behavior is consistent to the result in Ref. 25]. Regardless of the density-dependence of 



the symmetry energy, it is certain that the EOS of the neutron star matter is dramatically 
softened by the creation of the hyperons. Effect of the symmetry energy can be explored by 
calculating the mass and the radius of the neutron star, and comparing them with available 
observations. 

Figure [5] shows the mass-radius relation of the neutron star for np (Fig. [5(a)) and npH 
(Fig. [5(b)). In Fig.[5(a), NL3 shows the largest mass ~ 2.2M & which is well expected from its 
stiff est EOS. The difference from the other models is, however, not so much, approximately 
O.IMq. The maximum mass around 2M is a standard value which one can find in many 
works with the mean field theories. With octet baryons, the maximum mass drops to 
~ 1.6M (Fig. [5(b)). The difference in the maximum mass among the models is even smaller 
than that with the nucleons only, but the radius of the star at the maximum mass exhibits 
noticeable difference. For the FSU and the TW99 model, the radius at the maximum mass 
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is in the range 10.5 - 10.9 km, but the radius from the NL3 model is about 13.1 km. Precise 
measurement of the neutron star radius can give meaningful information and constraint to 
the nuclear symmetry energy at high densities. With the MDI model, on the other hand, 
we couldn't obtain a stable configuration of the star. It may imply that such a super- 
soft symmetry energy obtained from extrapolation from low densities is unrealistic at high 
densities, and one has to have more direct information from high-density region. 

IV. DISCUSSION 

In this work, we estimated the contribution of symmetry energy to the EOS of neutron 
stars with four different sets of parameterization. We also estimated the contribution of 
hyperons by employing DDRMF model without (np) and with hyperons (npH). In both 
np and npH cases, with soft symmetry energy, the neutron fraction increases as density 
increases at high densities, causing highly asymmetric dense nuclear matter. On the other 
hand, with hard symmetry energy, the neutron fraction decreases, and it makes the matter 
more symmetric. The dramatic change in the composition of particles in neutron star matter 
due to the density dependence of symmetry energy is very important to understand both 
the experimental results of heavy ion collisions and the properties of neutron stars. In other 
words, heavy ion experiments and neutron stars can give constraints on the symmetry energy 
and the EOS of dense matter. 

Recent analysis on the n~ /tt + ratio in FOPI data of GSI heavy ion experiments [l| raised 
questions on the behavior of symmetry energy at high densities. They showed that very soft 
symmetry energy is favored to explain the experimental result. However, since the symmetry 
energy becomes negative at high densities, which is contradictory to the naive expectation 
based on the QCD symmetry, it still requires more careful investigation. 

Other possibilities which we didn't consider in this work is the phase transition to exotic 
states such as boson condensation or transition to quark matter. Since kaon is the light- 
est boson with strange quark, and kaons will compete with hyperons in the /3-equilibrium 
condition, kaon condensation at high densities can be important. The work with kaon 
condensation is in preparation. 

We confirmed that the soft EOS reduces the maximum mass of neutron stars. However, 
with the super-soft EOS in the MDI model, stable neutron star cannot be constructed. This 
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indicates that the super-soft EOS may not be realistic as long as we accept the standard 
general relativity. Most accurate estimation on the mass of neutron stars has been done in 
double neutron star binaries, and the neutron star masses in these binaries are < 1.5M Q , 
which are consistent with rather soft EOS. However, recently the neutron star mass in a 
white dwarf-neutron star binary J1614— 2230 has been estimated to be (1.97 ± O.O4)M 
This observation is very important because the mass estimation is based on the detection 
of Shapiro delay, i.e., time delay of signals due to the gravity of companion star. It is well 
known that Shapiro delay is one of the key observation to have very good mass estimation. 
If this number is confirmed by other observations, it can rule out many soft EOS, including 
super-soft EOS. Hence the validity of soft symmetry energy that gives soft EOS can be 
checked by future investigations on the neutron star masses. 
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